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1 INTRODUCTION 



Cosmological weak lensing I B artelmann & Schneider! l200ll) has 
evolved to be a valuable tool in cosmology. Weak lensing sur- 
veys have contributed significantly to the determination of the 
dark mat ter power spectrum and to t he estimation of its ampli- 
tude o-g JSclmeideret"ai]|l998t Ivan Waerbeke et alj|200cl) by the 
measurement of cosmic shear and have enabled the reconstruc- 
tion of the dark matter distributi on in rich clusters of g alaxies (e.g. 
iKaiser & Sauiresll993l:rSeitz et alJl99rilMenard et alj2003l) . 

So far, only static matter distributions have been considered 
but from the solution to Maxwell's equations in the framework of 
general relativity it follows that gravitomagnetic potentials gener- 
ated by moving masses should alter the predictions for light deflec- 
tion ISchneider et alj fl992). While gravitomagnetic corrections to 
lensing are small, being of order v/c, where v is the velocity of 
the deflecting mass, they may contribute to the weak cosmologi- 
cal lensing: The cluster peculiar velocities following from a cos- 
mological Af-body simulation like the H ubble-volume simulation 
( Colbe re et alJl200d: Ijenkins et alJl200ll) are well described by a 
Gaussian distribution with zero mean and a standard deviation of 
cr v — 300 km/s, which is a fraction of 10~ 3 of the speed of light. 
Thus, relativistic effects influence the lensing signal appreciably in 
i, 1% of all clusters. In filaments iColberg et al. 2004) where mat- 
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ABSTRACT 

Dark matter currents in the large-scale structure give rise to gravitomagnetic terms in the met- 
ric, which affect the light propagation. Corrections to the weak lensing power spectrum due to 
these gravitomagnetic potentials are evaluated by perturbation theory. A connection between 
gravitomagnetic lensing and the integrated Sachs-Wolfe (iSW) effect is drawn, which can be 
described by a line-of-sight integration over the divergence of the gravitomagnetic vector po- 
tential. This allows the power spectrum of the iSW-effect to be derived within the framework 
of the same formalism as derived for gravitomagnetic lensing and reduces the iSW-effect 
to a second order lensing phenomenon. The three-dimensional power spectra are projected 
by means of a generalised Limber-equation to yield the angular power spectra. While gravito- 
magnetic corrections to the weak lensing spectrum are negligible at observationally accessible 
scales, the angular power spectrum of the iSW-effect should be detectable as a correction to 
the CMB spectrum up to multipoles of £ 100 with the Planck-satellite. 

Key words: gravitational lensing, cosmology: large-scale structure, cosmic microwave back- 
ground, methods: analytical 



ter is funneled towards the clusters, velocities are even higher: Infall 
velocities up to a few 10 3 km/s have been measured. 

The integr ated Sachs-Wolfe (iSW) effect, or Ree s-Sciama 
(RS) effect ISachs & Wolfell 19671 iRees & Sciamll 19681) arises if 
CMB photons encounter time-varying gravitational potentials on 
their passage from the last-scattering surface to the observer. When 
transversing time- varying potentials, the energy gains and losses a 
CMB photon experiences in entering and leaving potential wells do 
not cancel exactly. In this way, one expects a net blueshift of CMB 
photons in forming voids and a net redshift in matter-accreting clus- 
ters of galaxies. 

The iSW/RS effe ct has been studied theoretically in individ- 
ual objects 1 Martinez-Gonzalez et al. 1990) and can be used for the 
investigation of cluster mergers 1 Rubino-Martm et al. 2004). More 
importantly, it is sensitive to mapping the large-scal e structure as 
it highlights the sites of active structure formation jkaiset 



M 
1 ( 



Martinez-Gonzal ez & Sanzj|l990l: | Martinez-Gonzalez et al 



1982; 



1992, 



1994; Seliaki ll996l) . Furthermore, the iSW-effect may turn out to 
be a powerful probe for dark energy's influence on structure forma- 
tion 1 Crittenden & Turok 1996), when combined with other trac- 
ers of structure. A nu merical approach has been undertaken by 
iTuluie & LagumJ <1995 JbT) . who followed photons through a cos- 
mological n-body simulation and carried out the line-of-sight inte- 
gration numerically. 

The aim of this paper is to determine the corrections to the 
power spectra of weak lensing quantities caused by gravitomag- 
netic terms and to derive the iSW power spectrum, both by apply- 
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ing perturbation theory. In comparison to preceeding treatments by 
ISelialdll996l) and lCooravl 12002), the novel approach taken to de- 
termine the iSW power spectrum is by relating it to the gravitomag- 
netic terms in considered in lensing. Gravitomagnetic corrections 
to len sing have indeed been observed by iFomalont & Kopeikinl 
(2003) in imaging radio waves from a quasar on Jupiter, which is 
an outstanding archivement in VLBI astrometry. Gravitomagnetic 
corrections to lensing in the large-scale structure would only be de- 
tectable by their n-point statistics or by topological measures like 
Minkowski functionals, that would be especially sensitive to the 
effect's intrinsic non-Gaussianity. Concerning the iSW-effect, there 
are a quite a few reports on its detection in WMAP d ata in cross cor- 
relation with various populations of tracer objects jA fehordi et al] 

I 2004 | Fosalba et all2003tlBoughn & Crittenderl2004tlNolta et alj 
200 4lHirataetall2004 . but so far it has not been possible to de- 
rive values for single multipoles based on CMB data alone. 

The paper is structured as follows: After a compilation of key 
formulae and the derivation of Limber's equation for vector fields 
in Sect.|5| the power spectrum of weak gravitational lensing is con- 
sidered and the correction terms due to gravitomagnetic potentials 
are worked out by perturbation theory in Sect. Then, the iSW- 
effect is related to gravitomagnetic lensing and its power spectrum 
is subsequently derived in a perturbative approach in Sect.|4| The 
results are summarised in Sect. [5] 



2 KEY FORMULAE 

The assumed cosmological model is the standard ACDM cosmol- 
ogy, which has r ecently been su pported by observations of the 
WMAP satellite 1 ISpergel et all2 003). Parameter values have been 
chosen as Q M = 0.3, Cl A = 0.7, H = lOOAkms -1 Mpc~' with 
h = 0.7, n B = 0.04, n, = 1 and cr 8 = 0.9. 



2.1 Structure formation 

The cosmic density field p given in terms of the dimensionless den- 
sity perturbation 6 = (p - (p))/{p), where (p) is the average density 
of matter. The 2-point correlation properties of the overdensity field 
6 are described by the power spectrum P(k): 



(S(k)S'(k')) = (2nfS D (k - k')P(k), where 
6(k) = J d 3 x S(x) exp(ikx) 



(1) 

(2) 



is the Fourier transform of the overdensity field S. The normalisa- 
tion of the power spectrum P(k) is given by the parameter cr 8 , i.e. 
the variance of S on scales of R = 8 Mpc/7?: 



r * = iX dkk2w2(kR)p(k) - 



(3) 



Here, W is a window function of top-hat shape, the Fourier- 
transform of which is given by: 



3 3 
W(x) = — [sin(x) - xcos(x)] = -J\(x). 

V-' V 



(4) 



The shape of the power spectrum P(k) cc k" ■ T 2 (k) is well approx- 
imated by the transfer functions T(k) suggested by Bardeen et al. 
1 1986). They read in case of adiabatic initial conditions: 



1 http://map.gsfc.nasa.gov/ 



, ln(l+2.34?) r , „ 
n « )= 2.34, t 1+3J 



' + {\6.\qf + (5.46qf + (6.71qf 



The wave vecto r k is commonly divid ed by the shape parameter 
r introduced bv lEfstathiou et al l ^^2]) for CDM models and ex- 
tended to models with Q. t 1 by Sugiyama ( 1995): 



klMvc- l h 

q(k) = with r = Q. M h exp -L1 B 



1 + 



V2A 



(6) 



In linear structure formation, each Fourier-mode grows indepen- 
dently and at the same rate. The time dependence of the overdensity 
field <5 can be described by the growth function D(a): 

d'(a) 



6{a) = S D(a) with D(a) = a 



rf'(l) 



(7) 



T he shape of d' ( a) is w ell approximated by the formula suggested 
bv lCarrollet"ai]< 19921) : 



d'(a) = -Sl M (a) 



.( 



2.2 Dark matter currents 

The continuity equation p = -div(pu) requires the existence of 
large-scale coherent matter streams j = pv superimposed on the 
Hubble flow due to the formation of structure. In Fourier space, the 
relation between density and velocity reads in the Eulerian frame 
in linear approximation: 



v(k) = -iaH{a)m)-p5(k) = -idf(Q)^6(k). 



(9) 



The 1 /^-dependence causes cosmological velocities to come pre- 
dominantly from perturbations on larger scales in comparison to 
those that dominate the density field. H(a) = dln(a)/dr is Hubble's 
function. The function / describes the dependence of the equation 
of continuity on cosmic time and is a function of the mass density 
n <PeebleJl980tlLahav et alJl99ll) : 



din 6 dln£»(a) 



: Sl(a)° 



(10) 



d In a din a 

In analogy to eqn. 0, time evolution of of dark matter current ve- 
locities in the comoving frame is described by G(a), 

g'(a) 



G(a) 



with g' (a) = H(a)f(Q). 



(11) 



The theory of peculia r velocity field s is reviewed in detail in Dekel 
119941) and lStrauss & Willickl <1995t) . 

In general, the effects considered here are sensitive to density 
weigted velocities. The Fourier transform of vector fields q(x) = 
S(x)v(x) can be derived with the convolution theorem: 



q(k) 



(12) 



(13) 



2 J (2n) 3 

where the integrand has been symmetrised in eqn. 1 1 31 . 
2.3 Limber's equation for vector fields 

For the derivation of the angular power spectrum of the gravito- 
magnetic corrections to weak cosmological lensing or that of the 
iSW-effect, a variant of Limber's equation is necessary that is able 
to deal with projections of vector fields q(x) instead of scalar fields. 
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The derivation presented here is generalised from lVishniad ll987l) . 
Consider a vector field q(x) and its Fourier transform q(k): 

3 k 



q{x) 



C d 3 k 
7TTT q(k)exp(-ikx) 



(14) 



2.3.1 Components parallel to the line-of-sight C\{{) 

For the power spectrum Cl({) of the components of q« parallel to 
the line-of-sight, one obtains: 



Any effect k in question is assumed to measure a projection of q(x) 
on the line-of-sight, where e is a unit tangent vector on the photon 
geodesic. W(w) is a general weighing function dependent on the 
comoving distance w which describes its redshift dependence and 
is later to be replaced by e.g. the lensing efficiency function: 

K = f dwW(w)[e-q] = f dwW(w) [e ■ q(k)] f exp(-ikx)(15) 
Jo Jo J (2tt) 3 

The decomposition of the projected field k(0) into spherical har- 
monics Yc„,(S) is: 



c»(0 



K(ff) = J]Tj k ^ Y ^ *"» = I dS1 ^tJfi) with ( 16 > 

/_n ,„-_(■ *J4x 



2£+l Ue-\m\)\ 
Y(m{S) = \ —, — a 77 — ; — tt.P cm (cos 9) exp (imtf>) . 
An V (t + \m\)\ 



(17) 



In the random phase approximation, one obtains for the variance 
{\Kf m \ ) of k{6) in two directions e t and e 2 : 



(\Kt m \ l ) 



f dwi W(wi) f dfi, Y em {ei) 
dw 2 W(w 2 ) dCl 2 Y- m (e 2 ) 

J4ir 



(18) 



d 3 k 



(2n) 



- expi-ike^v^expiike^v^deiqik)] [e 2 q*(k)]). 



According to the cosmological principle, there is no preferred ori- 
entation, which allows to replace (\a em \ 2 ) with its average value over 
all m for a given I: 



(19) 



The vector field q{k) can be separated into components parallel and 
perpendicular to the line-of-sight e: 

9 = 9n + 9± w i m 1\\ = e (l ' e ) ant l 9± = 1 ~ Q\\ = e x (9 x e )- (20) 

For the projections e ■ q ± = and e x q^ = are valid. Eqn. 1191 is 
further simplified by choosing the coordinate system in a way that 
the z-coordinate is parallel to the wave vector, e z || k. Introducing 
spherical coordinates (9, t/>) and putting x = cos 9 on obtains: 



— ^— fdkl^f dwiWiwi) f dw 2 W{w 2 ) (24) 
(2tt) 2 J J Jo 

dxi exp(-faiwi) I dx 2 exp(ikx 2 w 2 ) 
[xiPeo(xi)x 2 P m (x 2 )] (qn(k, Wi)q\{k, w,)>. 

The dx\ - and dx 2 -integrations can be performed by taking advan- 
tage of the connection be tween Bessel functio ns and Legendre 
polynomials lWatsonll952llAbramowitz & Stegurl 19721) : 



Jc(z) 



dx P((x) exp(i'zjc), 



which can be can be generalised to give: 

r +i i d" 

J dxx"P e (x)exp(izx) = -— J t (z). 

Inserting formula I26i for n = 1 yields the final result: 



(25) 



(26) 



(2n) L J J J 



dw 2 W(w 2 ) 



dv^i 



J[(kwi) 



d\v 



■Jf(kw 2 ) 



<^ll(A:,wi)^(A:,W2)>. (27) 



2.3.2 Components perpendicular to the line-of-sight C^{1) 

After reducing the summation to \m\ = 1, the power spectrum C^{1) 
of the components of q ± perpendicular to the line-of-sight reads: 

W = 7=, vl, n f d " 2 f&ViWiWi) f'dw 2 W(w 2 ) 
{2nyi{( + 1) J Jo Jo 



X 



J+l 



d.vi exp(-toiwi) I dx 2 exp(fa 2 w 2 ) 



f-Jl -xjPnix^yjl - xlPaix^iqJk^w^q'Jk, w 2 )). 

The integration over the polar angles xi and x 2 is slightly more 
complicated than the previous case. Inserting the definition of the 
associated Legendre polynomials P[„, for m = 1 gives another fac- 
tor of Vl - x 2 : 

dPi(x) 



Am p 

P tm (x) = (-l) ra (l-.x 2 )? —r^- 1 -> Pn(x) 



dx" 1 dx 
The derivative of the Legendre polynomial can be replaced via 



(28) 



9,1 • - = xq\\ and q ± ■ - = Vl - x 2 exp(-i(f>)q ± 



(21) 



Furthermore, with exp(iA:ew) = exp(/A:xvv), the expression for the 
correlator is separated into: 

(q(k) ■ q'(k)) = 



With these simplifications, the integrals over the azimuthal angles 
0i and <fc in eqn. I19i can be carried out. Inserting the orfhonormal- 
ity relation J^" d(f> exp [i'(n - m)($>] = 2nS„ m reduces the summation 
over m to a single term, which is m = for the components parallel 
to the line-of-sight and \m\ = 1 for the components perpendicular to 
the line-of-sight. The final expression for the power spectrum C K (l) 
is now split into the two orthogonal projections: 

CM = Cl(0 + C^(e). (23) 



(i - x 2 )— p e (x) = e [P ( -M - xP e (x)] , 

dx 



(29) 



and the integration be carried out by inserting relation I26i . Then, 
the two Bessel functions can be combined by using the Bessel func- 
tion's derivative relation: 



XiX 2 (qn(k) ■ q;(k)) + ^1 - x 2 <T* ^-^(q^qlm- (22) dz 



- [z J t (z)\ = z / M ( Z ) - — = - 



which yields the formula: 



dz 



, Je(z) 



f dx Vl -x 2 P n (x)exp(;zx) = {(€+ 1)- 
J-l z 

This relation allows the final result to be written as: 



(2n 

Ji(kw\) Je(kw 2 ) 



w 2 



p 

Jo 

{q^k^^q'Jk, w 2 )>. 



(30) 



(31) 



(32) 
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3 GRAVITOMAGNETIC LENSING 
3.1 Definitions 

Light propagation through a slowly moving perfect fluid can be de- 
scribed by an effective refractive index n cff which follows from the 
post-Newtonian expansion of the Raychaudhuri-equati on to sec ond 
order for a weakly perturbed space-time (Schneider et al. 1992): 



2 4 
n cff = 1 - — <D + — A ■ 



(33) 



Here, O is the scalar potential and A are the gravitomagnetic vec- 
tor potentials, e denotes a unit tangent vector along the photon 
geodesic. In this approximation, the metric takes account of the 
matter density p and the matter current densities j = pv (i.e. terms 
of order u/c), but neglects the stresses 7y = pi/jVj + pSjj. The small- 
ness of these terms (being of order i> 2 /c 2 ) makes them unobserv- 
able, but they would be sensitive to the velocity tensor v{Uj , i.e. to 
shear flows, velocity dispersions and turbulence. 

In the near zone of a system of slowly moving bodies the re- 
tardation can be neglected; in this case the expressions for O and A 
are given as solutions to Laplace's equation: 

P(V) 



AO(r) = 47iGp(r) 
AA(r) = 4nGj(r) 



^ <D(r) = -G 



I s ' 



\r-r'\ 



(34) 



A(r) = -G f dV-^-. (35) 
J \r-r'\ 

The dark matter flux j is defined as the momentum density j = pv. 

An expression for de/dw, i.e. the change in propagation di- 
rection, follows from the variational principle 6 f ds = 0. s de- 
notes an affine parameter. The deflection angle a, being defined as 
a = emitiai _ e flnai can be obtained by integration: 



2 r AC 
a = — I di V x O j I dsexrotA. 



(36) 



The derivative perpendicular to the line-of-sight is defined via 
V ± C> = Vfl)-e(e-V<I)). The first contribution to a in eqn. 1361 corre- 
sponds to the attraction g towards the deflecting mass via g = -VO. 

The second term, however, is due to the gravitomagnetic fields 
generated by the matter current densities j. This contribution is 
related to the dragging of inertial frames which gives rise to the 
precession of orbiting spinning tops in the particular case of ro- 
tation of the field-generating body (Lense-Thirring precession, to 
be measured by Gravity Probe B 2 ). This formalis m has been ap- 
plied to various astrophysical systems, namely by Ibanez 1 1983) 
to grav itational light deflection of a rotating galaxy and by[Sereno 
(2003), who considered light deflection on rotating stars. Further- 
more, corrections to the deflection angle in gal actic mic rolensing 
due to moving lenses have been evaluated by He vrovskvH2004l) . 

3.2 Gravitomagnetic lensing by the large-scale structure 

Adopting the Born-approximation, which states that the gravita- 
tional light deflection is weak such that the integral in eqn. 1361 can 
be evaluated along a straight line instead of the photon geodesic it- 
self, it can be seen that gravitational lensing is insensitive to deriva- 
tives of the potentials along the line-of-sight. Working out the de- 
flection angles a and the tidal matrix = dajjdxj while neglect- 
ing derivatives along the line-of-sight yields formulae analogous to 
the case of static lensing, but with the gravitational potential O re- 
placed by O- ^A|| . Thus, the sources of gravitational light deflection 

2 http://www.gravityprobeb.com 



are the matter distribution 6 and the component of the matter flux 
j\\ parallel to the line-of-sight. The gravitational light deflection is 
stronger, if an object is moving towards the observer, because the 
photon stays in the interaction potential for a longer period of time, 
and vice versa. 

With the source term 8 + one obtains for the 

lensing converg ence k up to the comoving distance w 
( B artelman n & Schneidej2001f) : 



k(G, w) = 



3/m 



r^^M^LM, (37) 

C~ Jo //c(Wmax)a(w") \ C j 



where f K (w) = w, if spatial hypersurfaces are flat, which is the case 
for £l M + fl A = 1- The redshift distribution of lensed population of 
background sources such as faint blue galaxies is described by the 
distribution p(z)dz, being recast in comoving distance, Z(w)dw = 
p(z)dz. Then, the average influence Z(w) of the lever arms of the 
optical path for a given configuration of source and lens is given 
by: 



Z(w) = dwZ(w) — -— — — • 
J„ Jk(W) 



(38) 



In this work, we assume the generic distribution in redshift z for 
faint blue galaxies fc.f. lEllisil997i) . 



1 1/3 

p(z)dz = p z 2 exp(-/) with — = -T 
Po P 



(39) 



with mean redshift (z) = T(4//3)/T(3/j3) - 1.5 and most likely 
redshift z mllx = (2/f3)'^ ^ 1.21 for p = 3/2. For the average con- 
vergence k, the final result reads: 



m 



I dvv Z(w)k(6, w) 
Jo 



3-ff 2 £2o C Wm ^ ^. /k(w)/_ 2 



2c 2 Jo 



dw Z(w) ■ 



For Z(w), the phenomenological fitting formula 
1 



Z(w) - Zo exp 



1 - [log(w/w )] 



(40) 



(41) 



with Z = 1.441, b = 3.186 and w = 2314 Mpc//i is used, which 
yields excellent agreement with the properly evaluated function, as 
shown by Fig.Q The fitting formula alleviates the need of numer- 
ically carrying out the integration in eqn. 1381 when projecting the 
dark matter power spectrum. 



3.3 Perturbative treatement 

When considering gravitomagnetic corrections to gravitational 
lensing, the source term 6 of static lensing has to be replaced 
by <7n = (1 + -V||)(5. It should be emphasised, that the fluctu- 
ations in a weak lensing shear field are predominantly caused 
by modes in fc- space, that are propagating perpendicularly to 
the line-of-sight ( Blandford et aljll99ll) . Evaluating the correlator 
(q±(k, w^q'JJc, w 2 )> yields apart from the dominating 2-point term, 



(q±(k, wi )q±(k, W2)>2 P i = D{ Wl )D(w 2 ){6{k)S'{k)), 



(42) 



contributions of 3- and 4-point terms. The 2-point term stated in 
eqn. 1421 is of order unity and is the basis of the conventional the- 
ory of static gravitational lensing. In the perturbative treatment, the 
coupling of A:-modes in nonlinear structure growth is neglected, in- 
tegrations are implicitly taken to be restricted to quasi-linear scales. 
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comoving distance w [Mpc//i] 



Figure 1. The redshift weighting function Z(w) (c.f. eqn. 1381 . rendered as 
a dashed line), and the fitting formula 1411 (solid line) in comparison. 



3.3.1 3 -point term 

The 3-point term (q±(k, wi)q* ± (k, W2))3 P t consists of four contri- 
butions and is of order v/c compared to the 2-point term (c.f. 
eqn. 

(qAk,w l )q' ± (k,w 2 )h pt = ^ J ^ | (43) 

(S(—k, Wi)v±(p, w 2 )S(k - p,W2)) + (6(—k,Wi)v±(k - p, w 2 )6(p, W2)}+ 

(6(k,w 2 )v ± (—p,wi)6(p - k,w t )) + (S(k,w 2 )v ± (p - k,wi)6(—p,w\)} 

Here, the relations d*(k) = S(-k) and v'(k) = v(-k) were in- 
serted, which hold for real fields. By using this fact, the condi- 
tion Yiik, = is fulfilled which allows the 3-point correlators in 
eqn. J43I to be reduced to products of 2-point correlators by virtue 
of the formulae derived in Appendix IaI This yields four terms of 
the type (u ± S)(6S)/c and two contributions (f±<5) 2 /c 2 of second or- 
der. 

The correlation function can then be projected onto a plane 
perpendicular to wave vector k: The component of the velocity 
in the celestial plane is given by v ± = k x (v x k)/k 2 and hence 
u± — v sin — v Vl - x 2 , with x = cos 9, where 6 is the angle of 
separation between p and k. In doing this, the contributions of the 
type (v ± 6) 2 /c 2 vanish because they contain a multiplicative term 
(S(k)u(k)), which is a vector field collinear to k. The remaining 
terms can be combined to give: 



(In) (q ± (k,w l )q' ± (k,w 2 )hpt 



An 



D( Wl )D(w 2 ) [g'Ot) + g'(w 2 )] 



J p 2 dp Jdx Vl-x 2 1 P(\p\)P(\p -k\)M(j>,p-k) 



1 



+ P(\k\) 



M(k, -p)^r + M(k,p- k) : P ji2 P(\p - k\) 



\P\ 



\p-k\ 2 



\P\ \p-k? 
(44) 



was inserted. The function M(p,p') is defined as: 



7 pp \ p' p 7 \ pp' 



(46) 



It should be emphasised, that this 3-point correlator does 
not take account of the evolution of non-Gaussian features in 
the correlation function (8{k\ )5(k 2 )S{ k^)) a n d their influence on 
len sing determined by Ijain & Seliakl ll997l) : iBernardeaul i 19971) 
and Takada & Jain 1 2003a b), which strongly affects weak lensing 
quantities on small angular scales. 



3.3.2 4-point term 

The last contribution to the weak lensing power spectrum 
evoked by gravitomagnetic corrections is the 4-point term 
(q±(k, w^q'JJi, W2)>4 P i, which is of order v 2 jc 2 and thus strongly 
suppressed. The derivation of the term is easy prey: It can be 
don e in complete analogy to that of the Ostriker-Vishniac ef- 
fect jQstriker & VishniaJ 1 986t IVishniad 19 87). where any optical 
depth depending on redshift needs to be replaced by the appropri- 
ate weighting function (c.f. Sect. I4.6t and conversions from dark 
matter densities into baryonic densites are to be discarded. 

The derivation evolves cross-terms between the velocity and 
density fields, perhaps the most elega nt way of red ucing it to a sum 
of 2-point correlations i s given by Ma & Fry (2002), using a result 
from lMonin & Yagloml il965albl) : 



(2n?( qi {k)q'(k)) Apt = Fi q (\k\) : 



i riv rpL^k-p-p, 

c- J (2nf J (2n) 3 



p- p' 



p p> 



fr,fn P ^P\'> P ss(\P'\) + rT,T-uP6v(\p\)P6v(\p'\) 

\p'\ \p J \ \p'\ \p'l\ 



(47) 



where the irreducible 4-point correlation Ps v s u {k) has been ne- 
glected. 

Following lMa & 13 [2002), the projection to be carried out is 
(27r) 3 (g ± (A)<^(A))4 pl = 2Y Ji j6iejP'q q (\k\), where e, and ej are unit 
vectors along the lines-of-sight. The expression for Pq q (\k\) is given 
by eqn. J47t . In neglecting the irreducible 4-point term one obtains: 



(2ny(q ± (k)q' ± (k))A ! , t 



c 2 J (2n 



(48) 



(1 - x 2 )P 6S (\k - P \)P m ,(\p\) - U , - ^ PsvQk - p \)P 5v (\p\) 

\k- p\ 

Inserting the time-evolution of the density-velocity and velocity- 
velocity cross correlation terms, 

<v x (fc,wiK(fc,w 2 )} = g'(w^g'(w2)(v x (K)v\(k)), (49) 
{uAk, Wl )S'(k,w 2 )} = g\ Wl )D(w 2 ){vAk)8'(k)}, (50) 

yields the final result: 

(2n) i {q ± (k,wi)q\{k,w 2 ))^ t = 4D(wi)D(w 2 )g'(wi)g'(w 2 )x 



dxP(\k - p\)PQp\) 



k(l - x 2 )(k-2xp) 
k 2 - 2xkp + p 2 



(51) 



In the integrand of eqn. 1441 . the replacement \p - k\ 2 = k 2 - 2kpx + 3.4 Corrections to the power spectrum 



p can be inserted. Additionally, the time evolution of the velocity- 

(45) 



density cross correlation function, 

(v±(k,wi)6"(k,w 2 )} = g'(w 1 )D(w 2 )(v x (.k)6*(.k)), 



The three-dimensional power spectra (q ± (k)q' ± (k)) of the matter 
currents parallel to the line-of-sight is shown in Fig.|5|for the var- 
ious n-point contributions. Compared to the dominating 2-point 
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Figure 2. Three-dimensional power spectrum (q±(k)q* L (k)) including dark 
matter currents perpendicular to the line-of-sight, split up into the 2-point 
contribution (solid line), the 3-point contribution (dashed line) and the 4- 
point contribution (dash-dotted line). Additionally, the 4-point term of the 
currents parallel to the line-of-sight (q^(k)q*(k)) is drawn (dotted line). The 
power spectra are given for the present epoch, i.e. a = 1 and z = 0. 



term, the 3-point term is smaller by more than two orders of mag- 
nitude on small scales, but it becomes important on large spatial 
scales beyond 1 Gpc where it levels out. On these large scales, how- 
ever, limitations due to cosmic variance and insufficient sampling 
due to galactic foregrounds cast doubt on its detectability. The lev- 
eling on large spatial scales of the 3-point term is due to the fact, 
that for small k all powers in p in eqn. <44> add up to zero, which 
results in a flat curve for (^ ± (A:)^^(A)) 3pl . In comparison to the 3- 
point term, the 4-point term is smaller by another three orders of 
magnitude, but in shape it very much resembles the 2-point term 
and its influence on the weak lensing power spectrum is safely neg- 
ligible. 



3.5 Projected lensing power spectra 

The final expression for (q ± (k)q* ± (k)) can be projected by means 
of eqn. iM\ to yield the angular power spectrum of any lensing 
quantity, for example the convergence k. The distance weighting 
function to be employed can be read off from eqn. 1401 : 



W h (w) : 



3HqC1q f K (w) 



2c- 



Jdw' 



Z(w )- 



(52) 



By substituting y = kw, the distance weighting Wl(w) can be com- 
bined with the time evolution of the correlators to yield the func- 
tions 

n 2 

(53) 



<Pi(k) 



<Pe(k) 



J" dyW L J -^D(y)G(y) J 7yW^ (|) ^Z»(y)(54) 

v 1 2 

£dyW L {^j^-D(y)G(y) , (55) 

which carry out the projection of the 3-dimensional power spec- 
trum (q ± (k)q* ± (k)) to the angular power spectrum C K {€) by convo- 
lution: 



fe(k)4 P t ■ 




comoving wave vector k [Mpc/h] 1 

Figure 3. Contribution tpc(k) of the 2-point terms (thick lines), 3-point terms 
(medium lines) and 4-point terms (thin lines) to the angular power spectrum 
C K (t) of the weak lensing convergence k as a function of wave vector k, for 
f = 100 (solid line), C = 300 (dashed line), t = 1000 (dash-dotted line) and 
C = 3000 (dotted line). The curves are normalised to unity. 



(56) 



CM = + 1) fdk (q x (k)ql(k)) x Vt (k\ 

(2ny J 

where the associativity of the time-evolution enables the 3-fold in- 
tegration in eqn. I32i to be reduced to a 2-fold integration. 

The functions ipf(k) 2p i, ipt(k)^ vt and <pe(k)4 pt are shown in Fig.0 
Clearly, the fluctuations on a certain angular scale described by the 
angular power spectrum C{{) are dominated by spatial fluctuations 
with a certain wave vector k, which leads the peak of the function 
<Pt(k) to shift with increasing £, The projection kernels tpc(k) for the 
different n-point correlation functions show the common feature of 
rising fast at small k, but their decays at large k vary appreciably, be- 
cause the increasing influence of the time evolution of the velocities 
G(w) makes the functions to drop faster with increasing values of 
k. In this way, the observed asymptotic behaviour is ifif(k 2pt ) oc k~ 2 
for the 2-point projector, but the (pdkspt) an d fe(k4 P t) exhibit faster 
decays that are not described by a mere power law. 

The angular power spectrum of the weak lensing convergence 
C K {£) and its corrections due to gravitomagnetic terms is depicted 
in Fig. |4] Even at the largest angular scales considered here, the 
function (p{(k) peaks at values of k at which the corrections of the 3- 
point term are negligible. The detection of corrections to the weak 
lensing power spectrum due to gravitomagnetic terms requires the 
measurement of weak lensing shear on very large angular scales, 
which is beyond feasibility with current technology. On large an- 
gular scales, cosmic variance additionally limits the observability 
of gravitomagnetic lensing. 



4 INTEGRATED SACHS-WOLFE EFFECT 
4.1 Definitions 

The growth of structure imprints additional anisotropies on the cos- 
mic microwave background (CMB) by the time variation of the 
gravitational potentials along the propagation path of a CMB pho- 
ton. This effect is called the integrated Sachs-Wolfe (iSW) effect 
in the regime of linear structure formation l Sachs & Wolfej]^)67[ 
lHu&Surivarm]|l994 and Rees-Sciama effect I Rees & Sciama 
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Figure 4. Angular power spectrum C K (f) of the weak lensing convergence k 
and its correction due to gravitomagnetic potentials. The contributions from 
the 2-point term (solid line), the 3-point term (dashed line) and the 4-point 
term (dash-dotted line) are given separately. 



1 1 968t ISeli all 1 9961 : ICooravl2002t) if the density perturbations grow 
nonlinearly. The perturbations AT of the sky temperature T can be 
written as: 



~T~~~J w ~dt' 



(57) 



where <5® jdt is the derivative of the gravitational potentials with 
respect to conformal time t. Similar to gravitomagnetic lensing dis- 
cussed in Sect. 13.11 (c.f. eqns. <33l and GD), the effect is of the 
order 1/c 3 . 



4.2 Connection to the gravitomagnetic potentials 

Using the definition of ®(r) and the equation of continuity, p + div j, 
where j = pv is the matter current density, one obtains for the time 
derivative of O: 

tur.A—n r#,wL = G r*,Tj& 

dt 



t J \r-r'\ J \r-r'\ 

Integration by parts with respect to d 3 r' yields: 

r , ,V -j(r') f,, , ,1 

G dV ^_1 = -G dV/(r')-V 

J r-r' J JK r-r' 



(58) 



With the identity 

, 1 
V' 



1 



(59) 



(60) 



\r - r'\ \r — r'\ 

the derivative with respect to the primed coordinate can be replaced 
by a derivative with respect to the unprimed one. In the last steps, 
the V-operator can be drawn in front of the integral and the defini- 
tion of A (c.f. eqn. <35» be inserted to yield: 

G / dV '" (/) - v F^ = - v '(" G / dV F^) = - dlvA - (61) 

Thus, the final result reads: 

-di\.4 — t = f dwdivA. (62) 



at 

Eqn. I62t has a lucid interpretation: The formation of objects such 
as clusters of galaxies requires the matter fluxes / to converge and 



to accumulate matter (p > 0). Consequently, potential wells deepen 
and give rise to the iSW-effect in regions where divA does not van- 
ish. The iSW-effect measures the rate of change of a potential with 
respect to conformal time, or equivalently, the divergence of the 
vector potential A, which is proportional to the accretion rate in the 
Lagrangian frame. 



4.3 Putting the Sachs- Wolfe effect in a cosmological context 

In order to relate the statistical properties of the iSW temperature 
fluctuations r{0) to those of the matter streams j(r), the auxiliary 
vector field \{6) is introduced, which is defined as the negative gra- 
dient of the iSW temperature fluctuation t(0): 



(63) 



i.e. x(0) points along the steepest descent in temperature from hot 
to cold patches in an iSW field. Inserting eqn. 1621 into the defining 
equation for x(6) and converting the derivation with respect to the 
angular variable into a derivation with respect to the comoving 
variable r by using V e = /k(w)V, enables interchanging integration 
and differentiation: 



(64) 



X(0)= i J dw /*( vv ) V ( divA ) = Jr J dwf K (w)AA. 

Additionally, the replacement V (divA) = AA is inserted, which 
is valid if rot rotA = 0. This is fulfilled in vorticity-free velocity 
fields, rot/' = 0. In linear theory, initial vorticity perturbations are 
damped and after a sufficiently long time, the linear velocity field 
should be completely irrotational. Even in the regime of quasi- or 
nonlinear structure formation, Kelvin's circulation theorem forces 
the flow to remain irrotational and described by a velocity potential 
until dissipative processes on smallest scales give rise to vortical 
flows. 

Inserting Laplace's equation in the comoving frame, AA = 
4nGa 2 (p)(Sv) with the source term j = Sv, allows to replace 
Newton's constant G and the ambient mass density (p) by using 



Pcii = 3Hq/(8jtG), (p) = n p crlt and <p> = <p) /a 



c J J a 



3H 2 Q 



dw 



AW j_ 

a(w) c 



(65) 



The structural similarity of eqn. 1651 with the weak lensing conver- 
gence eqn. I40i is striking. 

Now, the 2-point correlation of the iSW temperature gradient 
field^(S) is related to the matter flux density j(r). For the derivation 
of the correlation function C T {£) of the temperature field r{6) itself, 
one rewrites eqn. I63i in Fourier space, yielding: 



x(0) 



r d 2 i 

J (2* 



>x(t) = iem 



(66) 



The expansion into Fourier modes rather than spherical harmonics 
is permissible, because r is expected to show fluctuations on small 
angular scales, so that r can be considered on a plane locally tan- 
gential to the celestial sphere. Squaring immediately gives: 



CM) : 



l 2 C T {l) 



t(£+\)C T (e), 



(67) 



where the last step is a valid approximation for small angular scales. 
The complementarity of gravitational lensing and the iSW-effect 
and the lensing counterparts of iSW quantities are illustrated in the 
flow chart Fig. [5] 
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Ps(k) < density S(r) 
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Figure 5. Flow chart with correspondent quantities in gravitational lensing 
(left column) and the integrated Sachs-Wolfe effect (right column). The- 
quantity analogous to the iSW temperature fluctuation t(6) in weak gravi- 
tational lensing is the deflection angle a(0). 



4.4 Perturbative treatment 

In the following we adopt the approximation that the rate of 
change of a potential is constant during the photon passage and 
that the accretion geometry does not change significantly. In this 
approximation, the component dA z /dz of divA is cancelled by 
the integration in eqn. 1621 and makes the iSW effect to mea- 
sure the components of divA perpendicular to the line-of-sight, i.e. 
t oc div ± A = dA v /dx + dA T /dy. Consequently, the matter fluxes 
q ± (x) = j ± (x)/c = 6(x)v ± (x)/c perpendicular to the line-of-sight 
primarily give rise to the iSW-effect. Accordingly, the fluctuations 
in the CMB due to the Rees-Sciama effect, being sensitive to the 
components of the matter flux perpendicular to the line-of-sight, 
are dominated by the components of £-modes parallel to the line of 
sight. 

Power spectra of the form {q\\(k)q?,(k)) have been considered 
by many authors in the derivation of t he Ostriker-Vishniac effect 
(e.g. lVishniacll987lljaffe & Kamionkowskilll998l) . In order to ob- 
tain the projection onto the line-of-sight, (2n) 3 {q^(k)q^(k))4 pt = 

Yitj jpj -r^jj- ^P^dArl), has to be carried out, which can be interpreted 

as the quadratic form k T Pk with a unit vector k and the matrix 
P = Pq q (compare eqn. l47t . The matrix P introducing the scalar 
product k T Pk is positive definite, due to the reality of the density 
and velocity fields. 



(27t) 3 <^n(A:) I ?|J(A:)>4p, 



4 r d 3 p 

?J (2^)3 



x 2 P 6S (\k - p\)P m ,(\p\) + 



(k - px)x 
\k- P \ 



(68) 



PsuQk - p\)PU\p\) 



The scalar product pk is again equal to pkx, where x is the co- 
sine of the angle of separation. Inserting the velocity-density and 
velocity-velocity cross correlation functions with their proper time 
evolution, 



(v\\(k, wi)5*(k,w 2 )) 
(v {l (k, wiK(£, w 2 )> 



= g'(w 1 )0(w 2 )<u ll (A)<5*(*)>, 
= g'(wi)g'(w 2 )(u ll (k)vl(k)}, 

yields the final result: 

(27T) 3 (q ll (k,wi)q ll (k,w 2 ))4 P i = D(w l )D(w2)g'{wi)g'(w 2 )x 



(69) 
(70) 



2tt 



M 



+ 1 . kx(kx - 2px 2 + p) 

dxPM-pmw ^_ 2xk p p+ / . 



(71) 



4.5 Power spectrum of dark matter currents 

The three-dimensional power spectrum (qn(k)q^(k)) of the mat- 
ter currents perpendicular to the line-of-sight is given in Fig. [5] 
Its amplitude is by a factor of 4 smaller than the power spec- 
trum (q ± (k)q* ± (k)), because the iSW-effect measures the streams 
Sv in contrast to gravitomagnetic lensing, where the source terms 
read (1 + 2v/c)S. Despite the fact that different projections are 
considered, the shape and asymptotic forms of (q\\(k)q^(k)) and 
(q ± (k)q* ± (k)) are very similar. 



4.6 integrated Sachs- Wolfe angular power spectrum 

Fig.Qshows the angular power spectra C T (C) of the iSW-effect t(6) 
and C x (€) of the iSW temperature gradient x(6) which have been 
obtained by applying the projection formula 1271 to the spectrum 
(q\\(k)qUk)) with the weighing function 



3//5Q0 f K ( w ) 
c 1 a(w) 



(72) 



which can be read off from eqn. J65L The redshift-weightings and 
the time-evolution of the density and velocity fields can be com- 
bined, which yields the function I73i after substituting y = kw, 



4pt 



£'d> lSW g) 



dJ,(y) 
dy 



D(y)G(y) 



(73) 



which mediates between the 3-dimensional power spectrum 
{q\\(k)qt(k)} and the angular power spectrum C T (£) by convolution: 



(74) 



CtW= (2^F / dk{ m (k)q\(k))x^ ( (k)- 

Again, the 3-fold integration in eqn. <27> is reduced to a 2-fold 
integration. The shape of the function if/ ( (k) is depicted in Fig. [6] 
for various values of C. In contrast to the function (pdk) used in 
the projection of the lensing power spectra, the function if/[(k) is 
symmetric about its peak, which is caused by the replacement of 
Jc(y)ly with the derivative dj f (y)/dy. The fast variability is again 
due to the strong influence of the velocity time evolution G(y). 

The angular power spectrum C T (() of the iSW temperature 
fluctuations r{6) along with the primary CMB fluctuations and the 
limiting PZancfc-sensitivity is depicted in Fig.0 The angular power 
spectrum has an amplitude of ^ 3 x lO -11 K 2 at small C and shows 
a slow variation with the multipole order L The amplitude agrees 
well with the result from Seliak 1 1996), but the decline of the power 
spectrum on large angular scales could not be confirmed, which 
is due to the fact that for large angles, the Bessel functions 7/(x) 
are a poor approximation to the Legendre polynomials Pc(x). The 
position of the peak in the projection kernel tftf(k) suggests that 
on the largest scales considered here, the angular spectrum C T {€) 
is dominated by fluctuations at the maximum of P(k) on scales at 
At 1 = 10 Mpc. With increasing multipole order C, the peak in if/((k) 
shifts only slowly towards higher values of k, which explains the 
small variation of C x (£) = 1(1 + l)C T (f). 

The channel averaged Planck-sensilivy is described by iKnoxl 
ll995tlTegmark& Efstathi o"vlT996T) : 



Cnoise(^) 



Ancr- 

N ■ 
1 *nix 



expfe^+l)], 



(75) 



where =2 5.03 x 10 7 is the number of pixels and 6 b the FWHM 
extension of the Planck-beam. For the average amplitude of the 
noise o" c( r per solid angle subtended by a single pixel we use the 
quadratic harmonic mean over all HF/-channels: 
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Figure 6. Contribution i//{(k) of the 4-point term to the angular power spec- 
tram C T (C) of the iSW temperature fluctuations r as a function of wave 
vector k, for t = 100 (solid line), t = 300 (dashed line), I = 1000 (dash- 
dotted line) and t = 3000 (dotted line). The curves have been normalised to 
a peak value of unity. 
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Figure 7. Angular power spectrum C^tV) = Tq MB C t (£) of the iSW tem- 
perature fluctuations r(0) (dashed line). The CMB power spectrum Ccmb(^) 
for the ACDM model (solid line) and the limiting P/anc/:-sensitivies 
CnoiseC) for angular resolutions AS = 5'0 (dash-dotted line) and A8 = 9' 1 
(dotted line) are depicted for comparison. 



n-2 Zj rr 2 



cr eff =13.42 /jK. 



(76) 



The sensitiviy considerations suggest that the iSW-effect is 
well above the noise level of the combined Planck HFI-channe\s, so 
that the power spectrum of C T (£) should be observable for angular 
scales 1 1 200 as a contribution to the primary CMB fluctuations 
Ccmr(I) , which i n Fig.0have been computed using the CMBfast 
code by Seliak & Zaldarriasa 1 1996). The signal can of course be 
amplified by cross-correlation with a suitable population of tracer 
galaxies. 



5 SUMMARY 

The scope of this paper is to derive the corrections to the power 
spectrum of weak gravitational lensing due to gravitomagnetic 
terms in the metric by perturbation theory. Within the same for- 
malism, the power spectrum of the iSW-effect can be determined 
as well. 

• The iSW-effect and gravitomagnetic lensing measure the evo- 
lution of velocities and densities in the large-scale structure and 
are sensitive to the cosmological parameters Q.m and erg. Applied 
to single objects like clusters, where the above described formal- 
ism equally applies, the iSW-effect would allow to measure the 
cosmological evolutio n of merger rates and dark matter accretion 
strengths Ivan den Bosch1l2002l : Iwechsler et af]l2002t Izhao et alJ 



1200: 

• Gravitomagnetic lensing would test general relativity on the 
largest scales (Mpc - Gpc) to second order in v/c, and could help 
to decide in favour of or against other metric theories of gravity. 
It should be emphasised that in the current theoretical description 
of structure formation or in current numerical simulations the mo- 
tion of bodies is described by classical mechanics, i.e. instanta- 
neous propagation of potentials and no relativistic increase of in- 
ertial mass with velocity, but the interaction of light with matter 
should be treated in the framework of the post-Newtonian limits of 
general relativity. 

• Gravitomagnetic terms influence the weak lensing power 
spectrum most notably on large spatial and angular scales, which 
are difficult to access experimentally. Furthermore, cosmic vari- 
ance and galactic foregrounds prevent accurate measurements on 
the scales in question, i.e. i Gpc/h and above. The small grav- 
itomagnetic corrections could be amplifie d by cross corre lation 
with the kinetic Sunyaev-Zel'dovich effect iSunvaev & Zel'dovich 
Il972l) . once future CMB telescopes will provide accurate measure- 
ments of line-of-sight velocities or with the velocity information 
from optical galaxy surveys. For contemporary weak lensing sur- 
veys, gravitomagnetic corrections to the cosmic shear do not play a 
significant role. 

• The iSW-effect is described by a line-of-sight integration over 
the divergence of the gravitomagnetic potentials. By this argument, 
the iSW-effect is reduced to a second order lensing effect. Every 
iSW quantity has a correspondence in weak gravitational lensing 
and the derivation of the power spectrum C T (C) proceeds in com- 
plete analogy to that of any weak lensing quantity, for instance 
that of the convergence C K {£). The most important difference of the 



derivation presented here to the ones carried out by Seliak 1 1996 
or ICooravlfcOoi) is that our derivation explicily pays tribute to the 
lensing nature of the iSW-effect. 

• Gravitomagnetic lensing and the iSW-effect are complemen- 
tary in measuring the matter flows parallel and perpendicular to 
the line-of-sight. The picture emerging is that (subject to the ap- 
proximations made) in gravitational light deflection (including the 
gravitomagnetic term A z ), the photon's A-vector is rotated but its 
normalisation is conserved. Contrarily, the components of A trans- 
verse to the line-of-sight change the normalisation of the A-vector, 
i.e. the photon's energy, but leave the direction of k invariant. 

• Both effects, gravitomagnetic lensing and the iSW-effect, are 
achromatic which makes them only accessible by their n-point 
statistics. Furthermore, the iSW-effect needs to be separated from 
other achromatic CMB structures such as the kinetic Sunyaev- 
Zel'dovich effect and the Ostriker-Vishniac effect. The derivation 
predicts iSW temperature fluctuations of a few jiK. on large angular 
scales, which is within reach of future CMB experiments like the 
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European P/aiic/c-mission 3 - 4 , provided that the modelling of Galac- 
tic foregrounds is sufficiently accurate to access these large angular 
scales. 

• The gradient ;^(#) of the iSW temperature fluctuation field r(6) 
should directly map regions of large matter flows, e.g. filaments 
and clusters with high peculiar velocities, but it can be expected to 
be very susceptible to noise due to the differentiation required in 
obtaining x(0) from t(0), which is reflected by the fact that ratio of 
the angular power spectra C x (f)/C T (f) is proportional to 1(1 + 1). 

The verification of the theoretical approach by a ray-tracing 
simulation of photons through a cosmological n-body simulation 
will be the subject of a future paper. The non-Gaussian features the 
iSW-effect and gravitational lensing exhibit and the mode-coupling 
in nonlinear structure growth are unaccessible to perturbation the- 
ory and are important on small scales. The novel approach to the 
iSW-effect presented here should allow a much improved preci- 
sion in the numerical treatment, because inaccuracies in interpolat- 
ing the scalar potential's time derivative <9®/<9f for each integration 
time step and in integrating a rapidly oscillating function inherent 
the direct approach fe.g. lTuluie & Lae una 1995a b) are alleviated. 
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APPENDIX A: DECOMPOSITION OF MIXED 3-POINT 
CORRELATORS OF DENSITY AND VELOCITY FIELDS 

In order to evaluate the 3-point correlation function 
(6(ki)v(k 2 )6(k 3 )) in perturbation theory, the density- and ve- 
locity fields are decomposed into linear terms <5 (l) , and small 
perturbations 5 (2 \ w (2) : 

5(k) = 6 {i \k) + 8 (2) (k) and v(k) = v m (k) + u {2 \k). (Al) 

As shown bv lFrvlJl984 . the second order density perturbation can 
be written as: 

= f f P^(2nf6 D (p + p'-k)M(p,q)6(p)6(p' 



I 



d 3 p 



(2*) 3 



-M(p,k-p)6 w (pW\k-p) 



(A2) 



with the function M(p,p') being defined as: 

%/) = y4[^^HKf. (A3) 

7 pp'\p' pj i\pp'j 

Clearly, the function M is symmetric, M(p,p') = M(p',p) and 
has the properties that M{-p,-p') = M(p,p') and M(-p,p') — 
M{p, -p')- For the first order perturbation of the velocity field, one 
obtains: 



v (2) (k) = -iH(a)m)-6 (2 \k). 



3 http: //planck.mpa-garching .mpg .de/ 

4 http: //astro. estec. esa.nl/Planck/ 



(A4) 



' http://www.cmbfast.org 



The 3-point correlation function (6(ki)v{k 2 )S{k 3 )) can now be ex- 
panded to yield to second order: 

<tf(*iM*2)*(*3)> - (S {l \kM l \k 2 )6 {2 \k 3 )) + (cycl) + 0(2) (A5) 

with the zeroth order term {S'- >> (ki)u l - [ \ki)S'- r> (ki)) vanishing due 
to v (l \k) oc (5 (1) (&) for truly Gaussian random fields. If the pertur- 
bation is contained in the density field S, inserting eqn. <A2t into 
the correlator yields: 

<^(fti)t/»(fe)d<»(t3)> = r x (^6) 

J (2ny J (2ny 

{2n?6 D {p + p' - k 3 )M(p,p')(S(k 1 )S(p))(6(p')v(k 2 )). 

Similarly, if the perturbation is the velocity-field v, one obtains: 

(S^WWKk,)) = f [pL x (A7) 

J (2?r) 3 J (2tt) 3 

{2n?6 D {p + p' - k 2 )M(p,p'){S(k x )v(p)){v(p')6{k 3 )). 

Collecting these results for the mixed 3-point correlator of density 
and velocity fields in question yields for the first order expansion 
of (6{k\)v{k2)S(ki)) in perturbation theory: 

(6 {1 \k,)v w (k 2 )6 {2) (k 3 )) = liikukdPisQkiWeAiJD, (A8) 
<<5 (1) (*,)!/ 2, (A:2)<5 <1, (A:3)> = M(k u k 3 )P Sv (\ki\)Ps v (\h\), (A9) 
(^(iiy^ft)) = M(k 2 ,k 3 )P S6 (\k 2 \)P Sv (\k 3 \), (A10) 

if the condition £f=i A, = is fulfilled. Hence, in first order pertur- 
bation theory, the 3-point correlation function can be decomposed 
into products of the density-density and density-velocity correla- 
tion functions, which are of the order v/c (eqns. lA8l and lA10l . and 
into the square of the density-velocity cross correlation, which is of 
order v 2 /c 2 (eqn. lA9t . 
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